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Abstract: 

The radiative corrections to the 3x3 scalar and the 2x2 pseudoscalar neutral Higgs 
boson mass matrices are calculated in the supersymmetric extension of the standard 
model including a gauge singlet superfield in the effective potential approach. The 
full t and b quark/squark contributions including the nonlogarithmic terms are taken 
into account, which are seen to affect the result significantly An analytic formula 
for the upper bound on the mass of the lightest neutral Higgs scalar including these 
corrections is given. 
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The supersymmetric extension of the standard model leads to constraints on 
the masses of the Higgs particles. Within the minimal extension the mass of the 
lightest neutral Higgs scalar is bounded from above by Mz at tree level. Recently 
radiative corrections to this upper bound have been computed by several authors 
[0J-0 with the result that contributions involving a large top quark Yukawa coupling 
can increase it by ~ 60 GeV. 

Within a supersymmetric extension involving a gauge singlet superfield []T0 1- |jT4 



the corresponding tree level upper bound depends on the value of a dimensionless 
coupling A similar to the non-supersymmetric standard model. Assuming, however, 
the absence of Landau-type singularities in the running coupling constants up to a 
GUT scale of ~ 10 16 GeV, A is bounded from above ||12||-||16|| implying a tree level 



upper bound on the lightest neutral Higgs mass of ~ 160 GeV. (The natural range 
of this mass for genuine choices of parameters has been estimated in fT/j| .) Also 
in this model radiative corrections to the upper bound involving the top Yukawa 
coupling have been computed using renormalization group techniques |I3]-P3 . 



These techniques, however, only take the logarithms of the ratio of the susy- 
breaking scale and the top quark mass into account, and since these logarithms are 



not very large, this is not necessarily a good approximation. Only in ||24|| , |p5|| also 
non-logarithmic contributions have been taken into consideration. 

It is the purpose of the present letter to compute the radiative corrections to the 
full 3x3 scalar and the 2x2 pseudoscalar mass matrices of the neutral Higgs bosons 
in this model using the effective potential approach. We include all contributions 
induced by the top or bottom Yukawa couplings, which are due to top/bottom 
squark/quark loops. (The other Yukawa and gauge couplings are genuinely smaller, 
and these contributions involve no colour factor N c = 3. At least in the minimal 
model they affect the final result only by a few GeV We end up with an 

analytic formula for the upper bound on the mass of the lightest neutral Higgs 
scalar, which includes the non-logarithmic radiative corrections. This allows a test 
of the leading log or renormalization group techniques, which fails significantly. A 



posteriori it justifies the procedure carried out in p5[ . (A similar formula for the 



upper bound has also been obtained by Comelli but there an expansion in the 
splitting between the two top squark masses has been performed.) 

The relevant part of the superpotential g of the model has the form 

g = h t Q- H 2 T C R + h b Q ■ E X B C R + \R X ■ H 2 S + ^-S 3 (1) 



where colour indices are suppressed and 



3 



Q = (b J ' Hl = {h[) ' H2 = {%) ' Q ' H = Ql ^ E ' etc - (2) 

The soft supersymmetry breaking trilinear couplings and masses are given by (sub- 
sequently the fields denote no longer chiral superfields as in (1), but complex scalars) 

(h t A b Q ■ H 2 T C + h b A b Q ■ H X B C + XAxH^S + ^A k S 3 ) + h.c. 
+m 2 1 \H 1 \ 2 + m\\El\ + m 2 s \S\ 2 + m 2 Q \Q\ 2 + m 2 T \T\ 2 + m 2 B \B\ 2 . (3) 
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We assume the Yukawa and trilinear couplings to be real. The scalar potential 
contains the standard F- and D-terms, the soft susy-breaking terms and in addition 
the one loop radiative correction of the form 

V rad = -^rStr[M 4 ln(M 2 /Q 2 )]. (4) 

(In the literature slightly different formulas can be found which are related to each 
other and the present one by a redefinition of the renormalization point Q 2 .) Next 
we assume vevs of the fields Hi, H 2 and S of the form 

with hi, h 2 and s real. The equations for extrema of the full scalar potential in these 
directions in field space read 

hi[m\ + X 2 (h 2 2 + s 2 ) + 9 -{h\ - h 2 2 )} + Xh 2 s(ks + A\) + l -V^ d = (6) 

h 2 [m 2 2 + X 2 (hl + s 2 ) + 9 -(h\ - hi)} + Xhis(ks + A\) + \v; h f = (7) 

s[m 2 + X 2 {h\ + hj) + 2k 2 s 2 + 2Xkh x h 2 + kA k s} + XA x h ± h 2 + ]^V™ d = (8) 

with g = \(g\ + g|). The elements of the neutral scalar 3x3 mass matrix M 2 read 
in the basis (Re H®, Re H 2 , ReS) after the elimination of ml, ml and m 2 s using eqs. 
(6-8): 

Mlu = ghl - Xs^(ks + A x ) + An, 
ni 

h 

Ml 22 = gh 2 2 - Xs-±(ks + A x ) + A 22 , 

M 2 33 = 4/c 2 s 2 + kA k s — XAx^—^- + A 33 , 

s 

K12 = (2X 2 -g)hih 2 + Xs(ks + A x )+Ai 2 , 
M 2 13 = 2X 2 shi + Xh 2 (ks + A x ) + Ai 3 , 

M 2 23 = 2X 2 sh 2 + Xhi(ks + A x ) + A 23 , (9) 



with 



1 



fOT W> ( 10 ) 
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A,-,- = -\/ mc! - — V rad (11) 

where no sum over i is implied in eq. (11), and the index % attached to V rad as well 
as Vi denote hi, h 2 , s for % = 1, 2, 3. For the neutral pseudoscalar 2x2 mass matrix 
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M 2 , which is left after the Goldstone boson has been eaten by the Z-boson, we just 
give the tree level result (again after the elimination of m 2 ,m 2 and Trig): 

( \ 



Mi 



_\(A±±h^h+M X(2ks-A x )(h 2 + h 2 ) 1 / 2 



X(2ks - A x )(h 2 + h\fl 2 -Xh x h 2 Uk + v) - 3kA k s 



(12) 



Next we have to evaluate the radiative correction V to the tree level potential 
according to eq. (f|). To start with, we only take the top squark- and quark- loops 
into account. Thus we need the 2x2 top squark mass matrix of the top squarks 
Tr,Tl and the top quark mass in an arbitrary background H^H^ and S (here we 
anticipate the computation of the corrections to the pseudoscalar mass matrix). In 
the basis (T£,T£) the top squark mass matrix is given by 



Kj\H°\ 2 + m 2 



h t (A t H° 2 * + XSHl 



h t (A t H$ + \S*H°*) 
and the top quark mass by 



h 2 \H°\ 2 + m 2 



V 



h 2 \H: 



0|2 



The eigenvalues of the top squark mass matrix read 



(13) 



(14) 



m 



T!,T 2 



with w 



h 2 t \H»\ 2 + -(m 2 T + m 2 Q )±w 



m 



Q 



rn 



2 \2 



+ K 2 t \A t Hl* + XSHI 



r0|2 



1/2 



(15) 



In order to calculate the corrections Ay to the scalar mass matrix, we have to replace 
H® by hi, H® by h 2 , and S by s in ( |14|) and (|15D , and to insert the corresponding 
masses into the formula (f|) for V ra d- Then the Ay are obtained from eqs. (|T0|), ([11]). 
The results are most easily expressed in terms of the following expressions X, Y and 
Z: 



X 
Y 



3h 2 X 
32^2 



h\h\ + \(m\ + m\ 



w 



3hf (Xhis + A t h 2 f 
32tt 2 w 2 
3hj h 2 (Xh lS + A t h 2 ) 

167T 2 W 



\n{m% 1 /m% 2 ) + \n{m 2 ^ 1 m 2 T2 /Q A ) + 1 



w 



4J 



(16) 



ln(m| Jm? J 



The corrections Ay to the scalar mass matrix are then found to be 



An 



AhS -X + X 2 s 2 Y 



hi 
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A 22 = -^X + A 2 Y + 2A t Z + ^Hm 2 Tl mymt) 

ll 2 o7T 

A 33 = ~^^X + X 2 h 2 Y 

s 

A12 = A t sX + XA t sY + XsZ 
A 13 = (2X3^ + A t h 2 )X + X 2 sh 1 Y 

A 23 = A t h x X + XA t h l Y + Xh 1 Z. (17) 

The minimal extension of the sypersymmetric standard model can actually be ob- 
tained as a special limit of the present model in the form 

k — 0, A — ► 0, s — > 00 with Xs = fx fixed, AAas = 777,3. (18) 

In this limit the results for the upper left 2x2 submatrix of M 2 can be compared 
with the ones of ref. |J, and we found complete agreement. 

Next we turn to the radiative corrections to the pseudoscalar mass matrix. First 
we have to expand the quark and squark masses (14) and (15) around the vevs (5) 
into imaginary directions in field space: 

H\ ^hx + ilm(H^), H% ^h 2 + Um(H%), S + ilm(S). (19) 

This expansion has to be inserted into the formula (ffl) for V rad , and V rad , in turn, has 
to be expanded up to second order in Im(Hf), Imi^H^) and Im(S). This allows us 
to read off directly the radiative corrections to the 3x3 pseudoscalar mass matrix in 
the basis Im(H®), Im(H®) and Im(S). The radiative corrections have to be added 
to the tree level part of the pseudoscalar mass matrix in this basis, which we have 
not shown above for simplicity. After the elimination of m 2 ,m 2 and m| with the 
help of eqs. (6-8) it turns out that the Goldstone boson is given by the same linear 
combination of pseudoscalar fields as at tree level. It is possible to take account of 
the complete radiative corrections to the remaining 2x2 pseudoscalar mass matrix 
using the following simple rule: Wherever A\ appears in the tree level mass matrix 
given by eq. (12), perform the following substitution: 

XA\ -> XA\ + A t X. (20) 

Actually this substitution describes nearly the complete X dependence of the 
radiative corrections Ay to the scalar mass matrix as well. This concludes the t 
quark and squark induced corrections to M 2 and M 2 . As noted in ||, the explicit 
Q 2 dependence in these formulas is cancelled by the implicit Q 2 dependence of 
running couplings; in the present approximation it is the running of m\, A\ and the 
running of the couplings within |-Fff 2 | 2 induced by wave function renormalization of 
the superfield H2, which have to be taken into account. 

The radiative corrections induced by b quark and squark loops are most easily 
obtained after observing the symmetry of the model under (h t , A t , m 2 -,, T L , T R , H 2 ) «-> 
(ht>, Af), m 2 B) Bl, Br, Hi). Accordingly it is convenient to define twiddled corrections 
A Tj to M 2 by the rule 

A- ■ = Ay (/it, A t , m 2 T , h x , h 2 -> h b , A b , m|, h 2 , h x ) (21) 
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where the indices i,j are related to i,j by an interchange of 1 and 2. Now the full 
radiative corrections to M 2 are given by eqs. (9) with A^- replaced by A„ + A^- 
everywhere, and the A^ can be read off eqs. ([[7]) after using the rule (|21|) . Also the 
5-quark- and squark-induced corrections to the pseudoscalar mass matrix are that 



simple; instead of (20) we have to perform the substitution 



XA X -> \A\ + A t X + A b X 



(22) 



where X is defined analogously to Ay in eq. ( 21j) . 

In all cases we are able to compare the results with the ones obtained within the 
minimal extension in || after the limit (0), and find agreement. This concludes the 
main results of the present paper. Of course, the radiative corrections to M 2 and 
M 2 are not of great importance, unless some knowledge of the tree level parameters 
of the model is assumed. Such knowledge can be obtained if assumptions on the 
parameters at a big scale as the GUT scale are made, and the parameters are scaled 
down to the electroweak scale with the help of the renormalization group equations 
IT ]. In particular an upper bound on the coupling A can be obtained this way 
which, in turn, implies an upper bound on the mass of the lightest neutral scalar 
Higgs particle [P^|-[J25|. Our results allow to find the non-logarithmic radiative 
corrections to this upper bound, and subsequently we will restrict ourselves to the 
t-quark and squark contributions again. 

First we note that any diagonal element of the 3x3 scalar mass matrix M 2 con- 
stitutes such an upper bound. This holds also for any diagonal element of T M 2 0, 
where O is an arbitrary orthogonal matrix. Let us thus study T M 2 with 



O 



( \ 

cos (3 — sin (3 

sin (3 cos [3 



V 



tan/? = — -. 
hi 



(23) 



The middle diagonal element of T M 2 reads 



\2 

M 2 Z [cos 2 2(3+ — sin 2 2(3] 
9 

+ sin 2 (3[A 2 t Y + 2A t Z + 



3h 2 m 2 



\n(m 2 Tl m 2 T Jml)] 



+ cos 2 (3 A Vy + sin 2(3 As [ A t Y + Z] . 



(24) 



The first line of (0) contains the tree level upper bound of ||13|| , and the loga- 
rithmic term in the second line is the only one taken care of within a leading log 
or renormalization group approach to this problem (usually under the assumption 
m 2 ^ ~ m^ 2 ~ m,Q ~ m 2 -, [fL7H- |123|l , and note that the quantity X with its explicit Q 2 
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dependence has completely disappeared from (|24j). After an expansion in — m\ 2 
one recovers the result of ref. |24| up to terms involving gauge couplings.) 



In order to exhibit the numerical significance of the different parts of the radiative 
corrections to the first line in (0), we plot the upper bound on the mass of the 
lightest neutral scalar Higgs particle of (^4j) versus the top quark mass m t in fig. 
1 for the following typical choice of parameters: A = .5, A t = 500 GeV, mg=(500 
GeV) 2 , rr4=(300 GeV) 2 , s = 500 GeV. For tan/5 we choose tan/5 = 1 (full line) 
and tan/3 = 5 (dash-dotted line). The dotted part of these curves denotes the region 
where the value of the top quark Yukawa coupling h t exceeds 1.1, and thus h t runs 
into a Landau-type singularity below a GUT scale of ~ 10 16 GeV [16]. The tree 



level upper bound, the first line of (]24]) , is given by 87 GeV for tan/3 = 1 and 90 
GeV for tan/3 = 5. Also shown in fig. 1 is the upper bound, if only the leading log 
radiative correction is taken into account, as a dashed line (for tan/3 = 5). 

We clearly see the numerical difference between the leading log and the full result, 
which now can be obtained easily for any choice of the parameters from eq. (P4|). 
Accordingly the results of fL8|-[p3| on the upper bound of the lightest Higgs scalar 



have to be taken with some care. Justified, on the other hand, is the procedure in 
24j , p5|| in order to obtain the corresponding upper bound. 
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Figure Caption 

Fig. 1: Upper bound on the mass of the lightest Higgs scalar according to eq. ( P4| ) 
for a choice of parameters as indicated in the text. Full line: tan/5 = 1 (tree 
level result: 87 GeV), dash-dotted line: tan/3 = 5 (tree level result: 90 GeV), 
dashed line: leading log result for tan/5 = 5. Dotted part: h t exceeds 1.1. 
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